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INTRODUCTION  AND  SUMMiXRY 

TliiK  interim  report  on  eontract  N0014-74-C-0273  summarizes  work 
on  two  t.isks  ('omiirisinK  the  third  year's  effort.  The  work  was 
aeeomp  1 ished  by  various  task  leaders  at  the  Environmental  Researcli 
Institute  of  Miclii^an  under  the  overall  guidance  of  Principal 
Investigator  Fred  Thomson.  The  performance  period  for  the  work 
reported  herein  was  January  1976  - December  1976. 

This  introduction  contains  a description  of  tasks  and  a discussion 
of  their  relevance  to  the  beach  reccmnaissance  problem,  a summary  of 
progress  to  date  on  each  task,  and  recommendations  for  future  work. 

The  introduction  is  followed  by  detailed  technical  discussions  of  work 
on  each  task. 

1 . I DESCRIPTION  OK  TASKS 

Effort  on  two  tasks  was  continued  in  the  third  year  of  the 
contract.  The  purpose  of  each  task  was  to  build  on  a proven  capability 
’ ' lor  mapping  features  pertinent  to  beach  reconnaissance  using  remote 

sensing  technology  and  to  develop  new  or  refined  information  extraction 
I apab i 1 i t ies.  These  capabilities  will  yield  new  or  better  beach 
reconnaissance  information  from  remote  sensing  data. 

1.1.1  BEACH  ENVIRONMENT  TASK 

The  Beach  Environment  Task  (leader,  Robert  Shuchman)  is 
built  from  the  proven  capability  of  mapping  silicate  minerals  using 
reststrahlen  techniques  developed  by  Vincent  and  Thomson  [1].  In 
previous  years,  a data  base  of  the  spectral  directional  reflectance 
of  fifty  diverse  silica  beach  samples  was  assembled,  some  reflectance 
data  in  the  J-IA  iim  spectral  range  was  obtained,  and  empirical 
analysis  of  the  spectra  resulted  in  several  algorithms  which  showed 
high  correlation  between  linear  sums  of  wavelength-averaged  reflectance 
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• iiui  lit  ri'l  1 iT  t aiHi-  iMl  ioK  .iiui  imi'ortaiU  lu'.ich  pliysiial  ])ropcr  t i es  siicti 
.IS  mn  i s 1 11  re  conti'iit,  sizi’,  and  degree  of  sorting.  Analysis  of  the 

J- 1 1 ,.m  d.it.i  revealed  t lu'  extinit  to  wliieh  compositional  differences 
coultl  hi  ilisi  inguished  t rom  changes  in  the  silicate  res  t s t rah  I I'n  bands 
in  the  8-10  ,.ni  r.inge. 

bill  the  measuri'raent  ol  more  reflectance  spectra  of  beach  sands  in 
order  to  verily  the  utility  .iiitl  generality  of  the  algorithms  developed 
bv  empirical  analysis  seemed  1)  too  costly,  2)  not  very  complete,  as 
risults  would  111'  dependent  on  the  spectra  in  the  data  base,  and  3)  not 
conducive  to  an  unde  r.s  I .md  i ng  of  the  physical  processes  which  c.ause 
the  rel  lection  of  light  from  beaches. 

Consequently,  radiative  transfer  models  of  Gwynn  Suits  (2]  were 
.idapted  (under  KRl.'i  internal  funding)  to  calculate  the  reflectance  of 
beach  sand.  Supportive  measurements  lU'cessary  to  obtain  the  data  needed 
tor  input  to  Suit's  modil  were  carried  out  under  contract  finding. 

.Model  calculations  of  tin-  reflectance  of  a carbonate  .and  a silicate 
be.ich  sand  wire  madi'  and  results  .are  reported  in  detail  in  seition  2 
.ind  summ.arized  in  section  1.2 

1.1.2  BOTTOM  KKATURKS  MAPPING  TASK 

The  Bottom  Fe.atures  M.apping  Task  (le.ader.  Dr.  David  I.yzeng.a) 
is  built  from  a proven  lap.ability  of  water  depth  mapping  using  the 
selective  penetration  of  -water  in  various  spectral  bands  sensed  by  a 
p.assive  mu  1 1 ispec  t r.al  scanner.  hast  year  we  defined  suitable  spi'ctral 
b.ands,  defined  the  depth  penetration  capability  of  the  modified  ratio 
algorithm  (MRA) , and  began  to  determine  which  types  of  bottom  features 
wi're  d i sc r iminab 1 e , assuming  non-scattering  waters.  Also,  a model  was 
developed  for  calculating  the  radiance  from  shallow  scattering  waters 
where  the  reflected  energy  from  the  bottom  constitutes  a major  portion 
of  the  observed  radiance. 

This  year,  the  radiative  transfer  model  was  exercised  to  calculate 
the  spectral  radiance  of  three  different  bottom  types,  in  three  different 
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ivpi's  of  water,  for  a number  of  different  wavelengths.  Then  the 
perlurmance  of  the  MRA  was  evaluated  by  drawing  decision  boundaries 
optimally  separating  the  three  bottom  types  at  a "training  depth",  and 
computing  the  probability  of  misclassif ication  at  other  depths. 

Results  are  presented  and  discussed  in  section  3,  and  summarized  in 
s ec  t 1 on  1.2. 
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1.2  SUM>L\RV  OK  RKSUI.TS 

111  the  Uf.'U'h  Knv  i runiTifnl  Task,  tlio  following  accomplishments  were 
made  this  year: 

1.  A model  was  developed  (from  a vegetation  canopy  reflectance 
model  of  C.  Suits)  to  calculate  the  reflectance  of  beach  samples, 
given  the  composition  of  the  sample,  some  optical  properties 

ot  the  constituents,  and  a knowledge  of  the  beach  moisture 
profile  (on  a millimeter  scale). 

2.  Because  the  fundamental  optical  properties  of  beach  constituents 
were  not  available,  a program  was  completed  to  measure  the 
optical  properties.  Spectral  ahsorbtlon  and  scattering 
coefficients  and  the  forward  scattering  fraction  were  measured 
for  16  common  beach- formi ng  minerals  over  the  range  0.4-2. 2 gm. 

r 

3.  To  provide  knowledge  of  tiie  millimeter  scale  moisture 
variations  in  beaches,  some  field  studies  were  initiated  after 
it  was  found  that  adequate  information  did  not  exist  in  the 

1 1 terature . 

4.  Preliminary  model  calculations  of  spectral  reflectance  show 
good  agreement  with  the  measured  values  of  reflectance  for  the 
samt'  beach  samples.  A more  thorough  evaluation  of  the  model 
calculations  will  be  performed  before  using  the  calculated 
radiance  values  for  mul tispectral  algorithm  development. 

5.  The  preliminary  model  calculations  show  the  importance  of  iron 
stain  and  opaque  minerals  in  controlling  the  reflectance  of 
predominantly  quartz  beaches.  This  occurs  because  quartz  is 

a very  transparent  material  and  the  dominant  effect  of  the 
quartz  is  to  scatter  the  radiation.  Because  of  the  scattering, 
radiation  has  a greater  chance  to  be  absorbed  by  stain  or 
opaque  particles. 
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In  thf  Underwater  Features  Task,  the  fol Inwing  conclusions  wore 
reached: 

1.  An  evaluation  of  the  modified  ratio  algorithm  (MRA)  was 
conducted  using  a radiative  transfer  model  previously  developed. 
Effects  of  changes  in  wavelengths  selected,  changing  water 
type,  and  changing  deep  water  offset  signal  on  the  probability 
of  mi  sc  1 ass i f ica t ion  of  vegetation,  sand,  and  mud  bottom  types 
were  assessed. 

Changes  in  water  type  degrade  the  performance  of  MRA,  by 
reducing  the  range  of  depths  over  which  adequate  performance 
can  be  obtained.  For  example,  using  an  algorithm  optimized 
for  Jerlov  Type  3 water  and  deep  water  offsets  from  Type  1 
water,  50%  or  better  average  probability  of  misclass i f icat ion 
can  be  obtained  out  to  6 m water  depth.  The  same  algorithm 
operating  in  Type  5 water  achieves  this  performance  only  out  to 
3.6  ra  depth. 

2.  Changes  in  deep  water  offset  signal  from  that  appropriate  for 
the  water  type  degrade  performance  at  shallow  water  depths  but 
may  actually  Improve  perfornumce  at  large  water  depths. 

3.  The  major  limitation  to  the  MRA  technique  is  that  two  wave- 
lengths must  be  chosen  such  that  the  water  attenuation  coefficient 
is  the  same  in  each  band.  This,  coupled  with  the  need  to 

select  wavelengths  which  penetrate  the  water,  leads  to  low 
target  to  target  contrast  and  consequently  sub-optimum  per- 
formance for  targets  of  Interest  (e.g.,  sand,  mud  and  vegetation). 
The  reflectance  ratios  of  these  nviterials  arc  quite  similar  in 
the  bands  selected.  For  example,  the  average  probability  of 
misclasslf ication  among  the  three  targets  selected,  at  zero 
water  depth  was  13.2,  16.6,  and  26.7%  for  wavelength  pairs 
selected  for  water  types  1»  3>  and  5 respectively.  This  is 
marginally  acceptable  perforimmce.  At  zero  water  depth,  this 
performance  is  only  attributable  to  the  selection  of  bands. 
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There  are  more  promising  algorithms,  not  involving  the  ratio 
technique,  and  not  subject  to  the  requirement  that  wavelengths 
must  be  chosen  where  water  attenuation  coefficients  are  equal. 
The  perfornuncc  of  these  algorithms  should  be  investigated. 


t r .pMr  Rl  T WVIl  l OW  UllU  t AH'  MA  ll  >RI(  S T «<F  ‘ >Nl  . r P«I'T  > MK.Mi'.AN 


2 

I’OTKNTIAI.  USK  OK  ONK  FUNDKl)  DKVKKOPMENTS  IN  OTHER  PROORAMS 

Altiuiusli  llu‘  jirfSfnt  Underwater  Features  Mapping  Task  and  Beach 
Environment  Task  are  m<iking  important  contributions  to  basic  remote 
sensing,  we  wish  to  nott'  tiiat  tiic  developments  which  we  are  making  are 
also  potentiaily  of  use  in  on-going  6.2  programs  of  tlie  U.S.  Navy  and 
in  programs  of  other  government  agencies.  Some  of  tiie  anticipated  or 
potential  applications  of  these  techniques  are  discussed  below. 

2.1  Underwater  Features  Mapping  Payoff 

We  are  aware  of  several  programs  in  coastal  ecology,  bathymetry 
surveying  and  in  Offshore  Continental  Shelf  assessment  where  we  feel 
that  the  techniques  we  are  developing  under  the  Underwater  Features 
Task  could  be  applied.  In  coastal  ecology,  a NOAA  program  to  map  benthic 
algal  communities  in  Alaska  could  benefit.  Present  surveys  (performed 
by  ERI.M  in  the  summer  of  1976)  were  made  at  low  tide,  necessitating 
careful  scheduling  around  tide  and  weather  conditions.  In  Alaska, 
the  available  flight  time  because  of  these  constraints  can  be  very 
small.  Being  able  to  map  algal  communities  while  submerged  will  offer 
greater  operational  flexibility.  Additionally,  the  techniques  would  be 
lu’lpful  in  assessing  the  situation  of  oyster  habitat  in  Florida,  which 
we  know  is  of  some  interest  to  EPA. 

One  of  the  most  important  uses  of  the  underwater  features  mapping 
capability  is  in  extending  bathymetry  calibration  points  (obtained  by 
pulsed  laser)  to  calibrate  passive  mu  1 1 ispec t ral  bathymetry  algorithms. 

Such  techniques  are  of  great  interest  to  the  U..S.  Navy  for  Tactical 
Beach  Reconnaissance.  ERIM  is  currently  engaged  in  a multiyear  program 
of  developing  a breadboard  Beach  Reconnaissance  sensor /processor  package, 
under  funding  from  the  Naval  Coastal  System  Eaboratory  (Contract  N61339-77-C-C059) . 
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i.asi  , with  Ihf  im-reasoil  pressure  to  develop  offshore  oil  reserves 
ot  the  r.S.,  the  Bureau  of  I, and  Management  lias  been  iliarged  with  Llie 
ri’spons  ib  i 1 i I V to  survey  offshore  (iontinental  Shelf  areas  to  assess  the 
env i ronmiuU a 1 impact  ot  offsiiore  oil  drilling,  pipeline  construction, 
aiul  shore  proci’ssing  facilities.  The  underwater  features  mapping 
algorithms,  when  perfected,  will  serve  as  a useful  tool  for  mapping 
nearsliore  features. 

2.2  Bi-ach  Knvironment  Payoff 

Ihe  algtiritlims  for  determining  important  beach  parameters  (compo- 
sition, grain  size,  degree  of  sorting)  from  remotely  sensed  data  will 
be  usetul  in  addition  to  providing  basic  understanding  of  beach  forma- 
tion and  dynamics.  Information  obtained  can  be  used  in  traf f icabil ity / 
mobility  models  developed  by  the  U.S.  Army  Corps  of  engineers  Waterways 
K.xperiment  Stations  (Vicksburg).  The  t raf  f icab  i 1 i ty/mob  i 1 i ty  models 
in  turn  provide  information  essential  to  the  Navy  and  Marine  Corps  for 
Beach  Reconnaissance  operations. 

1 wo  other  applications  of  the  beach  algorithms  might  be  in  the 
location  ot  heavy  minerals  on  beaches.  This  is  important  to  companies 
mining  heavy  minerals.  For  example,  a large  amount  of  titanium  is 
obtained  1 rom  beach  deposits  in  North  Carolina.  Certain  rare  earth 
minerals  are  also  present  in  these  deposits.  Also,  by  mapping  grain 
size  distributions  on  beaches,  some  important  information  might  be 
ob’ained  on  energy  inputs  to  beaches. 
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BKACfi  KNVIRONMKN’T  STUDY 


5.1  Al’l'KOACH 

Tlif  pre.sent  effort  began  in  1974  . During  Llie  fir.st  two  years, 

SiJ  beai'li  samples  from  a broad  range  of  env  i ronmisits  were  collerted. 

The  i)hysical  and  chemical  cliaracteristics  and  the  spectral  reflectance 
in  the  visible,  near  infrared,  and  thermal  infrared  spectral  regions 
were  measured.  Analyses  of  these  measurements  demonstrated  correlations 
hi'tween  spectra!  reflectance  features  and  certain  chemical  and  physical 
propi'rties,  and  stimulated  efforts  to  develop  a tlieoretical  model  to 
expl.iin  these  correlations.  Later  this  data  will  prove  valualile  as  a 
test  set  lor  verification  of  this  model. 

S[)ec  i 1 i c.il  ly  in  the  first  two  years  of  beach  environment  ONR  work 
[1],  we  have  observed: 

(1)  that  near  infrared  (0.8-1. 8 gm)  reflectance  of  beach  samples 
decri'ases  mono  t on  lea  11  y with  moisture  content  by  volume.  There 
.are  some  residual  effects  of  grain  size  such  that  if  the 
s.ample  universe  is  stratified  into  coarse  (grain  size  1 mm) 

.and  fine  (gr.ain  size  • 1 mm)  the  correlation  between  moisture 
conti'Ut  and  reflect.ance  improves.  Sample  composition, 

tor  typiial  quartz,  quartz-feldsp.ar,  .and  some  heavy  miner.al 
be.aihes,  apparently  plays  a minim.al  role  in  deti'rmining  ne.ir 
inlr.ared  retlectance. 

(2)  The  <)u.artz  res t s t rah  1 en  l e.atures  at  8-10  vim  are  reli.able 
indicators  of  quartz,  more  or  less  independently  of  moisture 
contiait  (from  0-24%  by  volume)  and  gr.ain  size  (from  0.33-2 
mm  mean  grain  size).  The  n..ignitude  of  the  reststrahlen 
feature,  measured  in  ref  lectaive , is  proportional  to  tlie 
percentage  of  quartz  in  .a  sample.  Concentrations  of  10-15% 
feldspar  and  kaollnlte  can  be  d isc  r i ml  n.ated  in  the  reststrahlen 
features  as  a broadening  of  the  reflect.ance  te.ature,  particularly 
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• It  I 111'  longer  (10-11  ,.ni)  wavelengths.  Smaller  (i-'j-") 

I'om  ent  r.it  1 ons  of  heavy  mini'rals  ran  also  be  discerned  in  the 
spectra  as  ;i  v.iriatiim  in  the  reflet  tance  in  the  10-12  ,.m 
region,  beyond  the  tpiart z rest st rah len  bands. 

(1)  While  computer  statistically  selected  single  band  and  band  ratio 
models  (bands  in  the  visible  and  near  infrared)  seem  to  predict 
grain  size  with  a high  degree  of  confidence,  there  did  not  appear 
to  be  a rational  physical  justification  of  the  bands  or  band 
ratios  selected. 

The  current  effort  was  directed  toward  explaining  the  correlations 
which  wi*  observed  in  the  .analysis  of  measured  reflect.ince  d.ata.  To 

• iciomplish  this,  we  di'Vi'loped  theoretic.al  modelling  techniques  to 
compute  the  spei  tr.al  refU'ctance  of  beach  sands,  given  the  spectral 
properties  of  the  constituents,  the  iDmposition  of  the  beach  satid,  and 
the  moisture  profile  of  the  beach  sand.  Initially,  we  felt  that 

a number  of  these  parameters  would  be  available  from  literature, 

.Hid  we  could  exercise  the  developed  models  to  verify  and  improve  the 
•algorithms  developed  in  previous  work.  However,  careful  liter.ature 
review  reve.aled  that  the  required  data  on  opt  ic.al  parameters  of  beach- 
lorming  miner.als  and  the  moisture  profile  of  he.aches  (on  a scale  of 
millimeters)  did  not  exist.  These  dat.a  wi’re  obtained  .as  part  of  this 
yi'.ar's  effort,  and  they  represent  a substantial  contribution  to  the 
knowledge  of  the  optical  properties  of  beach- formi ng  minerals  and  of  the 
moisture  profiles  of  beaches.  But  .as  a result  of  the  need  to  measure 
this  b.asic  i nf  Drm.at  ion , the  model  exercise  was  not  begun  until  late  in 
the  year,  and  that  work  Is  not  complete. 

The  CCANOPY  and  AQUACAN  (leaf  reflectance)  models  for  Dr.  C.  Suits 
|2|  were  adapted  to  calculate  reflectances  of  be.ach  sands.* 

* 

The  adaptation  of  the  Suits'  models  to  the  ONR  effort  was  accomplished 
under  ERIM  IR&D  funds.  Dr.  Suits  consultation  time  on  various  aspects  of 
the  adaption  of  his  model  were  also  covered  by  IR&D  funds. 
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riK'  .ulapu.l  Suits'  imult'ls  (callt-d  SANDRliK  ami  A()CASANDRK F)  will  hu 
usiul  to  ditirmiuc  tin.’  practical  limits  ot  remote  sensing  algoritlims  tor 
determining  plivsical  and  cliemical  properties  ol  beaclies.  Tlie  new  models 
in-ed  tile  I'oe  1 1 ic  ien  t s ot  absorption  and  scattering  and  tlie  forward  scat- 
tering 1 faction  lor  each  mineral  comprising  tiie  beach  sand  (i.e.,  quartz, 
leldspar,  kaolinite,  etc),  the  average  numlier  of  grains  (particles)  per 
given  volume  from  whicli  average  cross  sections  of  eacli  miner. il  type  can  be 
comiuited,  .ind  tlie  moisture  de|Uli  profile  of  tlie  beach  to  calculate  tiie 
ret  lectance  of  tile  beach.  KRIM  using  its  Cary  14  spectral  rel lectometer , 
measured  the  hemisplieric.il  transmittance  and  reflectance  of  e.icb  be.ii  h 
lorming  mineral  at  ;i  number  ot  prescribed  thicknesses  to  obtain  the  luu'ded 
coefficients  and  forward  scattering  fr.iction.  The  measurements  of  tlu'  1 unda- 
mental  optical  properties  of  common  be.ich  forming  minerals  in  the  0.4-d.2  ..m 
range  of  the  e lei  t romagnet  ic  spectrum  .ire  the  lirst  such  measurimiuit  s 
reported  to  be  made.  The  results  of  this  phase  of  the  work  are  currentlv 
being  prepared  lor  journal  publication,  and  .ire  further  discussed  in 
Sii  t i on  d . 2 . 

To  obtain  moisture  depth  profiles  of  ty]iical  biuiclu-s  on  .i  millimeter 
sc. lie,  I.RIM  initiated  a c, ireful  search  of  the  liter.itiire  .ind  consul  ti'd 
v.irious  experts  in  the  field  of  soil  moisture  as  to  the  ava  i 1 ,ib  i 1 i t y of  this 
i nt ormat ion.  Tike  the  needed  optical  properties  of  minerals,  the  moisture 
d.ita  are  the  millimeter  profiling  scale  was  not  available.  KRIM  developed  .i 
measuring  apparatus  to  obtain  the  moisture  data  lor  use  in  the  AQUACA.V  model. 
In  situ  field  measurements  (along  Fake  Huron)  and  laboratory  measurements  were 
nude  to  obtain  volume  of  w.iter  as  a function  of  depth  from  the  top  of  the 

sand  surface.  The  laboratory  measurements  utilized  beach  samples  that  were 
previously  collected,  dried  and  then  wetted  to  simulate  actual  beach 
conditions.  These  derived  moisture  profile  inform. ition  was  used  in  the 
AQUA.SA.VDKKK  model. 

A laboratory  appar.atus  was  also  developed  to  empirically  measure  the 
"skin  depth"  of  passive  mul tispectral  techniques.  The  empirically  obtained 
skin  depth  measurements  were  compared  to  one  preliminary  "skin  depth" 
result  obtained  using  the  SANDRKF  model.  The  results  were  within  + 1 mm 
in  10  mm  agreement.  Those  results  are  discussed  in  section  2.3. 
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}.2  MODLI.INC 

3.2.1  IN'IRODIICTIO.N 
rhi'  i a|)ahility  for  dolormining  important  physical  proport  ios 

ot  tuMi  h s.inds  by  roiiutt  o sousing  tochni(|uos  dopimds  upon  tho  intor- 
.Ktion  I't  radiatio'i  with  tho  ''onstituont  imitorials  ol  tho  hoatli.  Sucli 
inli’faction  is  cortainlv  compli'x  but,  novor  t ho  1 oss , must  follow  tiu' 
laws  ol  n.ituro.  Tho  purpose  of  making  a mathematical  rofloctanco  modol 
of  this  coraplox  i nl or.ic t ion  is  to  achiovo  insight  into  tho  relationship 
hotwoon  tlu-  romi'loly  received  signals  and  the  physical  properties  of  tho 
beach  sand  that  are  ol  interest. 

A mallu-mat  ical  mt'dc'l  of  a physical  phenomenon  is  tho  result  of 
incorporat ing  tho  mathematical  expressions  of  the  laws  of  nature  as  they 
apply  to  a I'omplox  c i ri' urns  t anco  so  that  tho  conclusions  drawm  from  tho 
a^;somblod  expressions  correspond  to  the  physical  results  of  an  experi- 
ment under  similar  c i rcums  t ances  . Such  a model  requires  e’xperiment  al 
V. ilid.it  ion  using  .)  I i>w  but  divi’rso  circumstances  in  order  to  prove  that 
the  ossiMU-o  of  till'  phononnuion  is  contained  in  the  logic  of  tho  model. 

Once  v.ilidation  h.is  been  achii-ved,  the  model  mav  T>e  used  as  a readilv 
.iv.i  i 1 ,ibl  e .ind  inexpensive  substitute  for  experiment  under  all 
c i rcums I ances  within  the  scope  ot  the  circumstances  of  the  val iihit ion 
oxpi- r i men  ts . In  addition,  the  logic.il  structure  ot  the  modol  provides 
the  insight  into  tho  signific.int  interaction  processes  so  that  more 
geni-ral  conclusions  may  be  drawn. 

3.2.2  RKKl.KCTANCK  MDDKi.  CONCKPTS 

The  most  elementary  model  of  sand  reflectance  is  the  simple 
pl.ine  mixtures  model.  The  model  employs  the  assumptions  that  all  sand 
particles  are  opaque  and  ate  randomly  mixed.  The  surface  of  the  sand 
layer  exposes  sand  particles  in  proportion  to  the  product  of  their 
mean  cross  sectional  area  and  the  concentration  of  the  particles  In 
the  sand  mixture.  Thus,  the  reflectance  spectrum  of  the  mixture  is 

predicted  to  be  the  area  weighted  average  of  the  reflectance  spectra  of  ^ 
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thi-  v.iriouH  mint’r.ils  exposed  .it  the  s.ind  surf.ice.  Multiple  scattering 
betWL'eii  particles  is  .assumed  to  be  negli>'ible  .and  oni-  surface  particle 
is  .issumed  not  ot  obscure  from  view  an  adjacent  surface  particle. 

Ibis  elementary  model  fails  to  achieve  ^ood  accuracy  bec.ause  the 
transmittance  of  particles  in  .a  finely  divided  state-  may  be  quite  Large 
.and  multiple  scattering  of  radiation  between  p, articles  should  be 
significant.  In  addition,  the  packing  of  gr.ains  in  l.ayers  produces 
p.irti.al  exposure  at  the  surface  so  that  line  of  sight  to  some 
p.irticles  will  depend  upon  the  dirc-ction  of  view.  A more  complex 
model  is  required  to  aci'ount  tor  these  effects. 

The  more  complex  model  which  is  used  in  this  work  employs  the 
identic.il  concepts  that  are  employed  by  the  directional  reflectance 
model  tor  veget.ative  canopies  (d).  It  mav  be  hard  to  visualize 
ot  f-hand  that  the  interaction  of  radiation  with  a veget.ative  canopy  is 
homomorphic  with  the  inter.action  of  r.adi.ation  with  s.ands  because 
veget.it  ic'n  .and  sand  hardly  appe.ar  the  same-  to  the  eye*.  Nevertheless, 
the  essence  of  the  reflection,  t ransmiss ion , and  multiple  scattering 
phenomena  is  the  same  except  for  the  optical  parameters  of  the  components 
th.at  are  involved  as  long  as  the  wavelength  is  much  sma  1 1 er  th.an  the  p.irticles. 

The  first  assumption  of  this  model  is  th.at  the  scattering 
components  are  distributed  more  or  less  uniformly  in  horizontal  layers 
where  the  mixture  of  component  types  may  be  different  in  the  v.arious 
layers.  In  .a  corn  field,  for  example,  the  tassi  Is  always  appear  at  the 
top,  green  healthy  mature  leaves  appear  in  .a  middle  layer  and  necrotic 
leaves  appear  usually  near  the  soil.  In  sands,  the  action  of  wind 
and  waves  are  likely  to  stratify  mineral  mixtures  vertically  and 
certainly  moisture  content  varies  with  a vertical  moisture  profile. 

The  division  In  layers  is  done  In  order  to  quantize  st.at  ist  i cal  ly  what 
may  otherwise  by  nearly  .a  continuous  distribution. 

The  second  assumption  is  that  the  radiation  field  m.iv  be  divided 
Into  two  types  of  radiant  flux  — specular  and  diffuse.  The  specular 
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thix  r(.'pri‘Si.'iit  s t lu-  r.uliation  arriviiiK  from  llio  .mMino  with  ro<- 1 i 1 i near 
pri'paKa  L ion  and  passes  thr<iu(^h  the  holes,  eraeks,  and  voids  of  tlie 
ensemble  I'f  ramiomly  paekeil  eomponents  without  deviation.  The  diffuse 
t 1 ux  is  di'rivi'd  t torn  the  speeul.ir  Mux  <tnd  is  tli.it  part  of  the  specular 
tlux  which  h.is  been  intercepted  bv  .1  sc.it  t e r i n>;  component  at  le.isl 
onci'  .ind  is  scattered  in  both  forw.ird  nnd  backw.ird  directions.  In  ci 
vi'p,etative  I'.inopy,  spe.ular  I lux  trecpienlly  reaches  llu-  soil  level 
.ind  .ii'pears  .is  sun  flecks  on  the  soil.  In  a s.uui,  speeul.ir  flux 
diminishes  exponen t i .1 1 1 y to  negligible  proportions  in  only  a few 
millimeters  depth.  Optic.-iliy,  the  sand  is  ini  initelv  deep.  However, 
the  dittuse  flux  is  dirived  from  the  specular  flux  in  the  identical 
manni'r.  Ihe  diffuse  flux  m.iy  piuietr.ite  much  deeper  into  the  sand  than 
can  speiular  flux. 

Ihe  third  model  assumption  is  that  the  manner  of  scattering  by  mineral 
particulates  can  be  adequately  represented  by  replacing  each  mineral  parti- 
culate with  a set  of  equivalent  l.ambertian  panels  which  have  the  s.ame 
spectral  transmittance  and  ref lec t.ince  as  does  the  component.  This 
assumption  defines  a simplified  scattering  phase  function  which  permits 
one  to  c.ilculate  ensemble  reflectances  in  closed  form.  The  form  of  the 
sc.iitering  ph.ise  function  becomes  significant  when  single  scattering 
is  the  dominant  phenomenon.  .Sc;tt  t er  i ng  by  widely  dispersed  aerosols 
in  the  atmosphere  (e.g.,  smoke  and  dust)  exhibits  detailed  phase  function 
effects.  However,  as  the  degree  of  multiple  scattering  increases,  the 
det.iiled  features  of  the  scattering  phase  function  are  no  longer  signifi- 
c.int.  In  both  sands  and  vegetative  canopies,  multiple  scattering  effects 
domintite  because  of  the  high  density  of  scattering  components. 

The  fourth  model  assumption  is  that  the  diffuse  flux  moves  generally 
vertically  upward  and  downward  with  a Lambertian  angular  distribution  as 
.1  first  approximation.  Tht  reflectance  of  the  ensemble  is  calculated 
using  the  method  of  self-consistent  field  with  the  specular  and 
approximate  diffuse  flux  as  the  illuminant  of  components.  The  ensemble 
reflectance  is  not  necessarily  Lambertian  but  is  only  approximately  so. 

Both  vegetation  and  sand  meet  this  approximate  criterion. 
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Bi.1  aiLsi'  o I llu'  homomorpli  ii  ro  1 a I ionsii  i p of  radiativf  i nCerai-.t  ions 
willi  sand  and  vogi.' La  I i vo  canopios,  tfu'  ri- 1 I (.t  t anre  model  previously 
develi'ped  for  ve^;eLaLion  is  applied  to  sand  with  the  appropriate 
eompoiu-nt  prtipert  ies  for  sand  inini-rais  suhstituted  for  veKet.itive 
e anop i es . 

3.2.3  Oi’TICAl.  PKOi’KRTlKS  OF  SAN'D  COMPONFNTS 

I’he  cross  section  view  of  a hypot  het  i cn  1 sand  layer 
containing  twt)  kinds  of  mineral  is  shown  in  Figure  1.  The  irri’giiJar 
shapes  of  the  sand  grains  result  in  some  more  or  less  random,  loose 
paiking  with  manv  voids.  These  grains  are  to  be  replaced  by  a number 
ol  e<)uiv'alent  Lambertian  panels  which  will  intercept  approximately  the 
same  amount  of  radiant  flux  as  do  the  actual  grains.  The  spectral 
reflei  tance  and  t r.ansm  i 1 1 ance  of  the  panels  are  to  be  the  same  as  the 
spectral  reflectance  and  transmittance  of  the  grains.  Since  these 
spectral  properties  may  cliange  with  grain  size  or  state  of  division, 
some  means  of  calculating  the  appropriate  spectral  properties  is 
retju  i red  . 

An  auxiliary  reflectance  and  t r.insmi t tance  model  for  mineral 
thin  sections  was  developed  in  order  to  relate  the  inherent  spectral 
properties  wliich  are  characteristic  of  a mineral  to  the  properties 
of  that  mineral  in  any  state  ot  division.  Tliree  characteristic  bulk 
properties  were  taken  to  he  sufficient  for  this  purpose  — the  spectral 
absorption  coefficient,  a,  the  forward  scattering  fraction,  FS , and 
the  scattering  coefficient,  s. 

Radiation  which  penetrates  a mineral  m.iy  be  absorbed  and  converted 
into  heat  energy  depending  upon  the  chemical  composition  of  the  mineral. 
A mineral  wliich  is  internally  homogeneous  without  inclusions  and  cracks 
will  transmit  radiation  passing  through  it  in  accordance  with  the 
relat ion. 
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Layer  1 


l^ayer  2 


Cross  Section  of  Sand  Layers 


Equilivant  Lambertian  Panels 


Layer  1 is  shown  as  consisting  of  grains  of  only  one  kind  of  mineral. 
l.ayer  2 is  shown  to  consist  of  a mixture  of  two  kinds  of  minerals. 

The  equivalent  Lambertian  scattering  panels  are  illustrated  on  the  right. 


FIGURE  1 . CROSS  SECTION  OF  SAND  LAYERS 
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wlirro  l^(x)i.s  tlu'  irr;uli;uK’c  on  a plane  at  depLii  x in  tlie  miner.il, 
ii  is  the  irradiaiu-e  on  a plane  inside  the  first  surface, 
a is  tile  spectral  absorption  coefficient. 


ilu'  spectral  absorption  coefficient  will  be  a function  of  the  wave- 
length ot  the  penetrating  radiation  and  will  depend  upon  the  chemical 
composition  of  the  mineral.  The  spectral  absorption  coefficient  is 
largi-ly  responsible  for  the  spectral  variations  in  mineral  reflectance 
.iixi  t !\insmi  ttanrc. 

The  forra.it  ion  of  miner.ils  is  a complex  natural  process  so  that 
miner.'ils  may  not  be  optically  homogeneous.  Foreign  materi;ils  are  often 
formed  in  the  interior  of  the  mineral.  Such  features  as  fractures, 
gas  bubbles,  small  crystals  of  associated  minerals,  and  grain  boundaries 
of  anisotropic  crystals  create  inhomugenei t ies  within  mineral  bodies. 
These  inhomogeneities  reflect  or  scatter  and  deviate  penetrating 
r.id  i .at  ion  from  rectilinear  prop.agat  ion.  A collimated  beam  of  radiation 
which  propagates  rect i I inearly  through  the  boily  of  a mineral  will  be 
diminislu'd  due  to  such  scattering  by  the  relation. 


K(x) 


K 

o 


(2) 


where  K(x)  is  the  irradiance  of  rectiline.ir  (lux  on  ;i  plane  at  depth  x, 

K is  the  irradiance  of  rectilinear  flux  .at  the  first  surface, 

o 

s is  the  scattering  coefficient. 


The  radl-ation  which  is  scattered  will  tend  tti  propagate  either 
deeper  into  the  miner.al  — forward  scatter  .and  contribute  to  trans- 
mi  tt.ance  or  reverse  direction  and  prop.agate  back  out  of  the  mineral  — 
back  sc.atter  .and  contribute  to  reflectance  of  the  mineral.  The 
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Ir.iction  ol  Si- a 1 1 o rt'd  radinLion  wliLch  coiRimifs  (k>i‘pc‘r  Into  t tio  mineral 
is  the  torward  scattering  fraction,  KS.  The  fraction  of  scatti'rc'd 
radiation  which  reverses  is  then  (1-FS).  Tlie  scattering  coefficient 
and  tile  torward  siattering  traction  can  be  spectrally  dependent  hut 
should  not  he  tlie  primary  determinant  of  the  spectral  cpiality  of  a 
mi nera 1 . 

ihesi'  three  optical  properties  of  a mineral,  a,  s,  and  FS , are 
assumed  to  be  independent  of  tlie  thickness  of  the  mineral.  That  is, 
the  Inhomogeneous  structure  and  chemical  composition  of  a mineral  thin 
section  is  assumed  to  be  evenly  distributed  so  that  a,  s,  and  FS , are 
inherent  properties  characteristic  of  the  kind  of  mineral  and  not  the 
size  of  the  mineral  sample. 

The  auxiliary  reflectance  and  transmittance  model  for  mineral  thin 
sections  makes  use  of  these  three  properties  and  the  index  of  refraction 
of  the  mineral  to  yield  the  thin  section  transmittance  and  reflectance 
for  any  mineral  thickness.  The  value  of  a,  s,  and  FS,  for  each 
wavelenj^th  must  be  determined  experimentally.  However,  these  properties 
cannot  be  determined  by  direct  experiment.  Instead,  a spectrometer 
is  used  to  measure  the  transmittance  and  reflectance  of  mineral  thin 
sections  havinj^  various  thicknesses.  Since  a,  s,  and  FS  are  presumably 
thickness  invariant,  the  value  of  a,  s,  and  FS  may  be  determined  by 

tindiiiK  the  value  ot  a,  s,  and  FS  which,  when  used  in  calculating 

tile  t ransmi t t ances  and  reflectances  of  thin  sections  having  these 
various  thicknesses,  yield  matching  values  for  the  reflectances  and 
t ransmi t t ances  found  experimentally.  A computer  Iteration  technique 
was  used  for  this  purpose. 

Then  spectral  values  of  a,  s,  and  FS  for  some  of  the  common  sand 
minerals  have  been  tabulated.  For  certain  opaque  minerals,  such  as 
hematite  and  limonite,  the  values  of  a,  s,  ind  FS  could  not  be 
determined.  'I'hc  transmittance  and  reflectance  of  these  minerals  are 
independent  of  grain  thickness  for  all  thicknesses  that  are  likely  to  be 

found  in  sands.  The  reflectance  and  transmittance  of  pure,  unstained 
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rliMi"  (|u,irtx.  is  liiu’  almost  oiiliroly  U)  siirtaoo  oftocts  wliicli  aro  also 
i nili'|)omU-iU  ol  Ufa  in  thioknoss. 

rlio  o|)tical  jiropi’ rt  i os  of  sand  rompononts  aro  dot  i' rmi  nod  and  aro 
introdnood  into  t ho  saiul  ro  t 1 oi- laiu’o  model.  llio  spectral  t ransm  i t t atii  e 
and  ri’tloctaneo  o’  the  I'qnivalont  Lambertian  panels  I’or  each  mineral 
are  determined  nsin>;  the  auxiliary  rel’lectanee  ami  transmittance  modi'l 
tor  mineral  thin  siations  where  the  thickness  ol  thin  section  is  the 
mean  urain  thickiii’ss  tor  each  mineral.  I'he  mean  cross  secticm  of 
ur.iins  ot  I'ach  miner.il  type  is  multiplii’d  by  the  cor  ri'spimd  i ng  number 
of  suih  gr.iins  per  u.iit  volume  for  a given  beach  sample  and  represents 
till’  siattering  el  feet  (' f the  equivalent  LambiTtian  panels. 

1.2.4  .modi; I.  SOl'TWARL  DKSCR  1 l‘T  1 ON 

I'igure  2 indicates  the  flow  of  information  beginning  with 
the  r.iw  experimental  CARY  transmittance  (t)  and  reflectance  (i  ) data 
and  ending  with  the  predicted  reflectance  spectrum  for  sand.  As 
indicated  in  soi  l ion  2.2.  3 the  neetled  input  into  the  SANDRKK  and  AQl'ASANDRKF 
modi'l  are  transmittance,  and  reflectance,  ( , v.ilues  for  the  indi- 
vidual miner, ils  in  .1  particul.ir  s.ind  conf  igur.it  ion.  These  t and 
v.ilues  must  be  c.iliulated  for  eai  h p.irticli’  size.  I'.xpe  r i men  t ,il  1 v , 

; and  . specti.i  were  determined  for  each  mineral  at  given  th  i cknessi-s . 

The  minerals  selected  are  shown  in  Table  1.  .Some  mineral.s  which  ware 
ojiaque  exhibited  a negligible  amount  ot  t ransmi  t t .nice . Iron  stains 
Were  treated  as  separate  miner.ils  .ind  i’X[)e  r i m«-n  t .i  I ' and  . spectra 
Were  determined  for  them.  These  tiles  ari’  repri-sented  in  the  tirst  row 
ot  Figure  2. 

In  order  to  predict  t and  , 's  f'M  any  p.irticle  size,  properties 
independent  of  size  can  be  determined.  These  projterties  previously 
discussed  on  the  forward  scattering  (FS),  scattering  coetficient  (s) 
and  the  spectral  absorption  coefficient  (a).  The  program  ABSCAT  (see 
Figure  2)  calculates  these  coetticients  for  a giver  mineral  using  the 
experimental  ’’s  and  t>'s,  the  experimental  thicknesses  and  the  parameters 
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COMMtlN  BI'.ACH  FORMING  MINFRAI.S 

(jiiarl  7. 

Foiy-crysta  1 1 int‘  (Quartz 
Milky  Quartz 
Ort  luH'  1 ase 
IMa^ioc  last* 

Kao  1 i n i to 
llomat  i t e 
Ch 1 or i t f 
01 i vene 
Miiscov  i to 
I'ornblomio 
Py roxono 
Garnet 
Magnet  i te 
I 1 men i te 
Carbonate 


FIGURE  2.  MODEL  FLOW 
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RI’A  itnii  RAP  (ri’lati'cl  to  the*  index  of  refraction).  The  experimental  data 

is  first  smoothed  to  eliminate  anomalies  due  to  the  experimental 

procedure.  The  equation  for  this  smoothing  is  as  follows  for  * to  > . 

1 n 


(■,)  Co)  r(*  ) 

'(■o)  = + -2-  ^ -4- 


c I _ 1^)  ( ■ ( ■ 14  1^)  O) 

___ — , 


( .J 


for  all  >j  = Lo 

The  next  program  Is  PREDICT.  This  program  uses  the  scattering  and 
absorption  coefficients  and  the  parameters  RPA  and  R^\P  to  predict 
t's  and  p's  for  a given  thickness.  These  predicted  values  are  now 
ready  to  be  used  to  model  the  sand. 

Theoretically  all  minerals  can  be  processed  through  the  ABSCAT 
and  PREDICT  programs.  But  as  mentioned  previously  some  minerals  cannot 
be  modelled  this  way.  Opaque  minerals  fall  this  model  because  the 
transml ttances  are  nearly  zero.  The  program  OPAQUE  smooths  the 
experimental  reflectances  and  averages  these  p's  for  all  thicknesses 
used.  The  transmlttances  are  set  to  zero  for  the  whole  spectrum.  An 
Important  assumption  made  is  that  these  average  reflectances  are  the  same 
for  all  sizes  of  particles. 
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St.iiiis  (f.K-t  iron  stain)  .iro  ilso  troatoci  as  minerals.  Tlie  stains 
oannot  tu-  proc ossisl  in  the  AKSCAT  ami  I’KKDICT  programs  because  these 
programs  rtsiuire  .it  least  two  expe  r i mint  t .a  1 tliicknesses  to  m.ike  .a 
predict  ion.  Tlu'  expi' r imen t .a  1 's  .ami  's  were  only  determined  .at  one 
thickness  lor  iron  st.ain.  Thus  the  program  STAIN  simply  smooths  the 
I'xper  imen  t a 1 v. lines.  The  .issumption  to  he  kept  in  mind  is  th.at  these 
smoothed  values  hold  for  all  thicknesses  of  st.ain. 

■file  progr.im  SNADRKK  com|nites  the  reflect.ance  of  mixtures  of 
miner. ils.  'flu-  inputs  to  this  program  consist  of  tlie  predicted  's 
ami  , 's  tor  tin'  minerals  constituting  the  sand,  the  amount  of  iron 
St. tin,  the  depth  ot  t hi-  s.and , the  reflectance  of  an  intinite  depth 
hai  kground  underneath  the  sand,  the  particle  cross  sections  .and  p.acking 
'.actors  tor  the  v.arious  miner.als,  the  angle  of  illumination  from  a 
point  source  illuminator  (or  optional  Lambert  inn  source)  and  the  angular 
position  o!  the  sensor. 

The  output  of  SANDRKF  is  a predicted  reflectance  spectra  of  the 
s.ind  tor  v.arious  cont  igurat  ions  of  illuminator  position  ( ) and  .sensor 
position  (!,•).  fhe  SANDRKF  program  was  derived  from  a program  which 
predicted  ri- 1 1 i-c t ances  for  a two  layered  vegetation  canopy  with  a back- 
ground. Thus  the  s.and  can  h.ave  two  layers  witli  different  der.sities 
and  grain  sizes  in  each  layer.  A background  material  for  the  sand  can 
also  be  specified.  fhesi-  properties  of  .SANDRKF  allow  great  flexibility 
for  modeling  v.arious  conf  igur.at  ions  such  as  thin  sands,  sand  after  a 
storm  with  a layer  of  organic  matter  or  stones  on  top,  etc. 

AQUAS/\NDREF  is  also  indicated  on  Figure  1.  Although  not  fully 
operational  at  the  present  time  this  model  is  very  similar  to  SANDRKF 
and  essentially  requires  the  same  inputs  with  the  .iddition  of  moisture 
information  on  a millimeter  scale.  Appendix  A of  this  report  includes 
a description  of  all  the  ONR  developed  beacii  environment  computer  programs. 
Included  with  the  description  of  the  program  are  the  experimental  derived 
i,  s,  and  FS  for  one  of  the  common  beach  forming  minerals  us*’d  in  thi 
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stmlv.  The  i,  s,  .uul  I'S  v.iliu-s  ol  thesf  mitUT.ils  .iri’  t lu‘  first  rt’pfirft'd 
v.il  Ill's  in  till'  I).  ..m  r.mni'. 

i.J.5  I’Kl.l  IMiNAKY  Mnl)i;  1. 1 N(i  RF.Sin.TS 

I'wi'  hisu  Ill'S  wiTi'  si'l  I'lti'ii  tor  pro  1 irai  nary  ovalnafion  of  tlio 
SANDKi  1 iiiiiiii'l,  oiii'  l arhonati'  and  t lii'  otlior  prim  ipally  qnarfz.  Tlu' 

I arlioni  1 1'  saiin'li'  '*as  lollciti'd  t rom  the  ocean  side  at  Marathon  Ki'v  , 
I'lofida.  Ihc  sample  is  comprised  of  rounded  and  plate  carbonate  shell 
tra>;ments,  calcareous  al^ae  flakes,  and  a trace  of  black  carbonate 
debris.  I'i^ure  1 shows  a reflectance  curve  in  the  0.4-2. 5 ;.m  region  of 
a s.impli'  oI  the  he.ich.  The  sample  was  measured  using  the  TRIM  Carv  14 
spe  trori'f lectometer .The  curve  of  the  actual  measured  reflectance 
ippiars  ,i.s  a solid  line.  Superimposed  on  Figure  3 appears  the  generated 
SANDRFF  rel  li'ctauie  of  the  c.arbonate  beach.  This  curve  is  dotted. 

Ihe  results  .iri'  Very  encouarging,  for  although  the  absolute  value  of  the 
lurves  ditfer,  the  overall  shape  and  absorption  bands  seem  to  correspond 
remark. ibly  well.  After  comparing  the  gener.ated  carbonate  reflectance 
to  the  CARY  ret  lect.ince  of  the  actual  sample,  the  origin.'il  beach  samples 
rel  I ect  .nice  in  the  visible  region  was  measured  bv  use  of  .a  photo-detector 
and  calibrateii  T.ambertian  rel  lectance  panel.  The  results  indie. ite  that 
the  Cary  value  tor  the  c.'irbonate  beach  s.imple  is  high.  An  ex.tmi  nat  i on 
ot  the  C.iry  d.it.i  is  currently  underway. 

A fpiartz  be.ich  was  also  used  in  the  preliminary  evaluations  of 
SANDKFF.  The  qu.irtz  be.ich  selected  w.as  t rom  Mexico  beach,  Florida 
(previously  di'Siribi'd  (Ilf.  Original  mineralogy  inform.it  ion  on  this 
be.ich  indii.ited  it  was  1007  quartz  with  no  iron  stain.  Wlieii  SANDRF.K  w.is 
run  under  this  .assumption,  the  results  proved  unsat  isf  ac  t orv  . Conse- 
quently, the  Mexico  be.ich  sample  was  reexamined  and  it  was  found  that 
slight  iron  st.aining,  black  organic  matter  (1%),  .ind  some  carbon. ite 
m.iterial  existed  in  the  sample.  SANDKFF  was  again  run,  this  time  with 
95.57  quartz,  2.5%  c.arbonate  1%  illmenite  (our  substitute  for  black 
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('r^;iinii  m.ittfr),  aiul  1%  iron  stain.  liu*  rfsults  arc  shown  in 
Ki>>iirc  A.  As  in  Fi^nrc  3 tiic  soliti  line  curve  rci>rcscnts  tiie  Cary 
obtained  rellectance  ot  tiie  actual  sample  while  the  dotted  line  is  the 
nenerati-d  SANDKKK  reflectance.  These  results  appear  in  close  agree- 
men  t . 
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WAVELENGTH  (MTCRONS) 


FIGURE  4.  COMPARISON  OF  MEASURED  QUARTZ  REFLECTANCE  WITH  SANDREF 

MODEL  RESULTS. 
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K i ASSKSSMK.N’T  OK  RyXOlAI'ION  PKNKTRATIUN  IN  BKACll  SANDS 
i.1.1  RATIONAI.K 

'nif  cloptii  to  wli  ioli  radiation  ponotratcs  beacli  samples 
(skin  di'pth),  and  depth  below  wbicli  chan^'es  in  beacli  parameters  do  not 
contribute  to  llie  radiance  viewed  by  the  sensor, are  of  utmost  importance 
imderst and i ng  the-  limitations  of  remote  sensing  techniques  for  beach 
reconnaissance.  First,  we  must  understand  which  portion  of  tlie  beach 
we  are  sensing  with  the  passive  techniques.  If  we  are  sensing  properties 
of  a layer  of  beach  less  extensive  than  that  needed  for  t ra f f i cab i 1 i ty 
analyses,  approaches  to  correlate  the  remotely  sensed  information  with 
the  traf 1 icabi 1 ity  modc'l  requirements  must  be  devised.  Second,  we 
wish  to  be  sure  that  the  2-1  cm  thicknesses  of  sand  that  were  measured 
in  tlie  Cary  14  spect  roref  1 ectometer  were  optically  tliick,  so  that  the 
ettect  of  the  Cary  14  sample  holder  reflectance  is  minimal.  Alternatively, 
it  the  samples  were  not  optically  thick,  we  would  have  to  derive 
appropriate  corrections  for  the  library  spectr.il  data  previously 
measured . 

».1.2  DESCRIPTION  OF  SKIN  DEPTH  MEASUREMENTS 

The  purpose  of  the  skin  depth  experiment  was  to  obtain 
some  estimates  of  the  m.jximum  depth  which  0.6-1.  1 pm  radiation  would 
penetrate  a sand  layer  and  be  back-scattered  sufficiently  to  affect 
the  r.idiance  seen  by  a passive  sensor,  and  to  verify  estimates  of  *^his 
depth  calculated  by  the  G.  Suits  model.  Results  of  the  G.  Suits 
SANDREF  model  indicate  that  the  skin  depth  for  quartz  beaches  should 
not  exceed  1 cm.  These  calculations  also  indicate  that  coarser  grained 
sands  and  sands  with  increasing  moisture  content  (above  15%  bv  volume) 
have  larger  skin  depths.  To  confirm  these  theoretical  calculations, 

• It  was  decided  to  obtain  empirical  measurements. 

An  experiment  was  desinged  to  mimic  both  the  CARY  14  optical 
configuration  and  potential  remote  sensing  geometry,  as  shi>wn  in  Figure  5. 
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riif  insiiic  ot  tlio  s.impK'  disli  was  lined  with  aluminum  toil  to  contain 
all  incident  li^hL  within  the  sample  and  avi>id  losses  through  the  side 
ol  till'  dish.  The  source  used  was  a GK  1323  tun>;sten  filament  densito- 
meter hulli  litted  with  an  aluminum  foil  cone.  The  sensor  we  used  was 
a Coiierint  Kadiation  Model  212  power  meter  couiiled  witli  a silicon 
detector  having  a 1.1  ,m  cut-off. 

file  lietector  was  fitted  with  a . fi4  3 ,.m  (deep-red)  C2-5B  Corning 
lilter.  Ihis  n-duced  the  spectral  band  of  operation  to  approximately 
.t)4 3-1.1  . m.  A flat  black  cylinder  was  fitted  to  the  detector  head 

to  provide  greater  d i rec t iona 1 i ty . A f ter  each  sand  sample  was  placed  in 
the  dish,  a llat  black  cardboard  ring  was  placed  on  top  of  the  sand  to 
minimize  extraneous  light  reflection. 

The  actual  measuremi'iit  of  skin  depth  was  achieved  by  sliding  a 
black  card,  then  a white  card  under  a petri  dish  filled  with  various 
thicknesses  of  sand.  If  fluxuations  of  greater  than  17  of  the  total 
incident  jiower  were  noted  on  the  power  meter  it  was  concluded  that 
light  was  penetrating  the  sand  layer.  More  sand  was  then  added  to  the 
petri  disli  and  the  test  was  repeated.  Wlicn  the  power  fluxuation 
was  reduced  to  1%  of  the  incident  (lower,  the  thickness  of  sand  was 
ri'Corded . 

Five  s.ind  sangiles  with  various  percentage  of  water  by  volume 
were  tested.  Samples  were  selected  which  v.iried  greatly  in  botli  grain 
size  and  mineralogical  composition  so  that  tiie  results  could  be  generalized 
to  all  the  sand  types  included  in  the  fifty  sample  spectral  reflectance 
data  base. 

3.3.3  RF.SULTS  OF  SKIN  DKPTH  MKASUREMHNTS 

Table  2 summarizes  the  skin  depth  measurements.  The  depths 
are  given  in  millimeters  as  a function  of  percent  moisture  by  volume. 

The  sample  numbers  are  the  same  as  those  in  the  previous  interim 
report  (l ) . 
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RKSULTS  OF  EMPIRICAL  SKIN  DEPTH  MEASUREMENTS 
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Till'  ri'siilLs  .in-  in  i-omp  1 ot  i>  aKri't-nuMU  with  Suit's  theorotical 
la  leu  1 at  ions  usinp,  the  moiiil  ied  SANDRK.K  Model.  The  measurements  show 
no  sipnitieant  penetration  beyond  8 mm  lor  all  the  sand  samples 
whi'thi'r  wet  or  drv. 

Skin  depth  imreases  with  water  eontext  as  predicted  by  Suits. 

Ihe  eo.irser  s.md  samples  KA I (2.67  mm  mean  grain  size),  CAl  (.81  mean 
gr.iin  size)  and  81)1  (.71  mm  me.an  grain  size)  h.ive  larger  skin  depths 
th.in  the  1 i in-  grain  samples  MX2  (.21  mean  grain  size)  and  HAl  (.82  mm 
me.in  grain  size),  again  in  agrement  with  the  theory.  The  empirical 
me.isurements  also  indicate  skin  depth  decreases  with  increasing 
pi'rcentage  ot  opacpie  minerals  and  of  iron  stain. 

1.4  IM.ANS  FDR  KinURF.  BFACll  F.NV  1 RON'MKNT  TASK  WORK 

As  discussed  previously,  the  collection  and  reduction  of  the 
mini'ral,  moisture,  and  .skin  depth  data  w;is  quite  rigorous  and  consumed 
a large  fraction  of  the  effort  allotted  to  this  task.  The  actual 
operation  and  verification  of  the  SANDREF  and  AtlUASANDREF  wore  not  begun 
until  the  final  quarter  of  the  current  year's  effort. 

For  next  year,  the  model  verification  begun  this  year  will  be  more 
vigorously  pursued.  The  verification  progr.im  will  be  completed  by 
calculating  the  reflectances  of  some  of  the  multi-mineral  sand  samples 
using  measured  physical  par.imeters  as  model  inputs  and  comparing 
these  with  the  CARY  14  measured  values  now  stored  in  the  data  bank. 
Further  software  refinements  of  the  SANDREF  and  AQUASANDREF  models, 
to  improve  efiiciency  of  operation,  will  be  made.  This  effort  will 
not  be  a large  one,  however. 

During  consistency  analyses  of  the  reflectance  and  transmittance 
values  of  the  mineral  samples  obtained  from  the  CARY  14  data,  we  noted 
that  the  reflectance  and  transmittance  values  totalled  greater  than  1. 
This  is  physically  Impos.sible  and  indicates  discrepancies  in  the  basic 
CARY  14  measurements.  The  problem  was  first  noted  in  the  measurements 
, for  iron  stain  and  was  subsequently  found  in  other  data.  We  feel  that 
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tliL'  nu'lhod  iisoci  Lo  moasurc-  llic  Lranstnission  of  samples  may  be  faulty  — 
some  radiation  incident  on  the  sample  nuay  be  leaking  around  the  sample 
and  impinging  directly  on  the  detector,  causing  transmission  values  to 
be  too  large.  The  cause  of  this  effect  will  be  more  thoroughly  explored 
next  year,  and  the  s.imple  data  corrected  as  appropriate. 

The  preliminary  exercise  of  the  SA.NDRKF  model  this  year  indicates 
the  importance  of  accurate  minera 1 og ical  information  about  the  sand 
samples.  We  feel  that  the  accuracy  requ  i rei.1  exceeds  the  accuracy  of 
the  data  obtained  on  our  fifty  sample  library.  To  provide  an  adequate 
model  evaluation,  we  will  repeat  the  mineralogical  analysis  for  certain 
of  the  samples.  Special  care  will  be  taken  to  account  for  trace  and 
opaque  minerals  in  predominantly  cjuartz  samples,  for  In  these  samples 
trace  minerals  appear  to  influence  the  spectral  reflectance  most 
dramat ically . 

Beach  samples  that  will  be  used  for  the  model  verification  will 
include  the  Atlantic,  Delaware  Bay,  Oregon  "heavy  mineral  beaches", 
Mexico  Beach,  Lake  Michigan  quartz-feldspar,  Kaolinite,  and  Marathon 
Key  carbonate  samples. 

Development  will  also  continue  on  improved  moisture  and  grain 
size  algorithms.  This  will  be  accomplished  using  the  proven  SANDRF.K 
and  AQUASA.MDRKF  m(jdels  to  calculate  the  reflectance  of  a number 
of  beach  samples  of  cliosen  composition  and  moisture  profile.  The 
composition  of  those  beach  samples  will  be  selected  in  such  a way  as  to 
verify  the  form  of  the  moistu’-e  algorithm  already  postulated  and  to 
attempt  to  derive  a more  meanin,  ful  mul tispectral  grain  size  estimation 
algorithm.  The  imporoved  moisture  algorithm  will  be  generated  in  such 
a way  as  to  be  minimally  sensitive  to  mineral  composition  and  grain 
size.  The  limits  of  performance,  accuracy,  and  skin  depth  ns  well  as 
the  optimum  spectral  channels  to  use  will  be  determined  for  both  the 
moisture  and  grain  size  algorithms. 
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UNDKHWATF-R  KEATURKS  STUDY 

rill.'  ob  ivc‘  fur  tills  year's  effort  was  two-told:  first,  to 

develop  a model  for  calculating  radiances  over  shallow  water  as  would 
he  measured  by  an  airborne  sensor,  and  second,  to  use  these  calculations 
to  evaluate  and  improve  data  processing  techniques  tor  mapping  bottom 
leatures.  Wt’  have  chosen  the  modeling  approach  because  it  allows  the 
simulation  of  any  set  of  conditions  and  is  thus  cciuivalent  to  collecting 
an  aitual  data  set  under  ideal  circumstances  where  all  the  relevant 
parameters  are  known.  Furthermore,  the  model  can  be  used  to  generate 
a data  set  with  no  noise,  representing  an  ideal  sensor,  or  with  an 
arbitrary  amount  of  noise,  simulating  the  data  collected  by  an  actual 
sensor . 

4.1  MODEI.  DEVKUJFMKNT 

in  order  to  calculate'  radiances  occtirring  in  the  atmosphere  above 
.1  body  of  water,  .a  water  radiance  model  was  developed  which  is  actually 
a combinatii'n  of  twc>  radiative  transfer  approximations.  In  this  model, 
the  double-delta  approximation  developed  by  Turner  111  at  ERIM  is  used 
to  calculate  solar  irradiances  and  atmospheric  effects,  and  an  exten- 
sion of  the  quasi-single-scatter  ing  approxi  m.i  t ii>n  developed  during 
last  year's  ONR  contract  (I)  is  used  to  calculate  the  water  reflectance. 
The  combined  model  includes  the  I'ffects  ol  atmospheric  sc.attering  of 
incident  sunlight,  reflection  of  light  .it  the  water  surface,  scattering 
.and  .absorption  of  light  by  the  water  I'olumn,  reflection  at  the  bottom, 
and  attenu.it  ion  of  the  reflected  liglit  in  tlie  atmosphere. 

A computer  program  was  written  in  FORTRAN  to  implement  this  water/ 
.atmospht re  model  on  the  University  of  Michigan's  Amdahl  470  computer. 

The  program  reads  in  four  sets  ol  input  parameters,  described  below, 

.and  computes  a set  of  radiances  as  a function  of  water  depth.  These 
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rad i anres , which  may  be  viewed  as  a simulaled  remote  sensing  data 
set,  are  stored  tor  later  display  and  analysis. 

ilie  tirst  set  ol  input  parameters  read  in  are  called  scene  para- 
meters, bocavise  they  describe  tlie  extemal  conditions  for  wiiich  the 
simulation  is  to  be  made.  Ttiese  include  tlie  atmospheric  visibility 

(or  meteorological  r.ange),  tlie  scone  elevation  above  sea  level,  the  solar 
zenith  angle,  the  atmospheric  scattering  functions,  and  the  wavelengths  at 
wliich  the  calculations  are  to  be  made.  The  program  then  calculates  the 
optical  tliickness  of  the  atmosphere,  based  on  Elterman's  [4]  tables,  and 
the  direct  and  diffuse  solar  irradiance  at  the  water  surface. 

The  second  set  of  input  parameters  describe  the  platform  from  which 
tile  radiances  are  observed,  including  the  altitude  of  the  platform  and 
tlie  view  angles.  After  reading  these  parameters,  tlie  program  calculates 
the  sky  radi.ince  and  the  patii  radiance  scattered  into  the  sensor  from 
the  atmospliere. 

Tlie  tiiird  set  of  input  parameters  are  tlie  water  parameters:  these 
are  the  diffuse  attenuation  coefficient  of  the  water,  the  volume 
scattering  function  of  the  water,  and  the  external  water  surface 
reflectivity.  The  latter  quantity  was  included  as  an  input  parameter  to 
allow  for  the  effects  of  variable  sea  state,  using  tlie  relationships 
established  by  Austin  [51.  The  surface-reflected  radiance  is  then 
computed  by  multiplying  this  parameter  by  the  sky  radiance  at  the 
appropriate  zenith  angle. 

The  final  set  of  parameters  read  include  the  bottom  reflectivity 
and  the  range  of  water  depths  to  be  considered.  These  are  read  last 
to  allow  duplicate  sets  of  calculation  to  be  made  for  different  bottom 
types.  The  water  radiance  is  calculated  using  equations  similar  to 
those  described  In  Appendix  (1  of  reference  1,  with  the  following  modifi- 
cations. First,  the  diffuse,  or  irradlanie,  attenuation  coefficient  (K) 
is  used  Instead  ol  the  beam  extinction  coefficient  (c)  to  describe  the 
water  attenuation.  The  approximate  relationship: 
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K = ) 

t'st.ib  1 istu’d  by  Cordon  ot  nl  (6  I is  nssiimod  to  exist  between  ttiese 
eoel  I Uients,  K Eieing  the  more  easily  measurable  and  more  widely 
reported  quantity.  Second,  the  contribution  from  internal  reflection 
at  tl)e  water  surface  is  computed  to  all  orders  by  assuming  for  this 
purpose  that  the  upward  radiance  just  below  the  water  surface  is 
isotropically  distributed.  This  is  a somewhat  more  accurate  approxi- 
ra.ation  th.an  the  first-order  correction  described  in  reference  1,  and  in 
particular  yields  the  correct  result  in  the  limit  of  zero  water  depth 
l/l,  for  a l.ambertian  bottom.  Third,  the  water  reflectance  is  computed 
separately  for  the  direct  (sun)  .and  the  diffuse  (sky)  components  of  the 
incident  irradiance.  The  sky  radiation  is  .assumed  t(5  be  isotropically 
distributed  for  this  c a 1 cul .at  i on . 

Raili.ance  calculations  .are  made  for  e.ai  h depth  specified  in  the  last 
set  ot  (i.ar.ameters,  and  stored  on  .an  output  device.  These  results  are 
displayed  and  an.alyzed  as  described  in  the  following  section. 

4.2  PARAMT.TKR  SKLECTION  AND  MODEL  RESUl.TS 

The  utiliz.ation  of  the  radiance  model  for  simulation  purposes  is 
presently  limited  by  the  .availability  of  b.asic  water  and  bottom 
reflect.ance  measirement s.  .lerlov  [8|  has  tabulated  spectral  irradi.ance 
attenuation  d.ata  for  various  water  types,  but  no  corresponding  scattering 
data  is  presented.  On  the  other  hand,  I’etzold  (9)  has  tabulated  complete 
sets  of  parameters,  including  volume  scattering  functions,  for  water 
samples  taken  from  a variety  of  locations,  but  only  at  one  wavelength. 

In  order  to  calculate  the  volume  scattering  functions  at  other 
wavelengths,  the  total  scattering  function  was  considered  to  be  m.ade 
vip  of  two  components:  a Rayleigh  component,  which  varies  inversely 

with  the  fourth  power  of  the  wavelength,  and  a particulate  component, 
which  was  assumed  to  be  independent  of  wavelength.  The  total  Rayleigh 
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scatloriiiH  coolticiont  .iL  0.55  tim  was  taken  to  be  0.002  m from 
(’eti'(jlci' s measurements  on  filtered  sea  water.  fhc  Kayleigb  component 
at  0.55  ,m  was  subtracted  from  the  total  scattering  function  at  this 
wavelength  in  order  to  obt.tin  tile  particulate  component.  Tiie  Rayleigh 
component  was  then  recalculated  at  eacti  wavelength  using  the  fourth- 
power  law  and  added  to  tlie  fixed  [larticuiate  component  to  yield  tlie 
total  volume  sctittering  1 uni  t ion  at  eacli  wavelength. 

•ierlov's  water  types  1,  1,  and  5 were  selected  as  representative 

of  the  conditions  typically  encountered  in  remote  sensing  in  the  coastal 
zone.  The  corresponding  scattering  data  sets  were  chosen  by  computing 
tlie  irradiance  attenuation  l oef  1 ic  ients  from  equation  (1)  for  each  of 
i’etzold's  water  samples  and  selecting  those  which  were  most  nearly 
equal  to  .Ierlov's  values  at  this  wavelength.  The  data  sets  from  Station 
9 (AUTEC),  Station  5 (HAOCK) , and  Station  2200  (NUC)  were  selected  on  this 
basis  as  corresponding  most  closely  to  water  types  1,  3,  and  5, 
respectively.  The  irradiance  attenuation  coefficients  for  these  three 
water  types  are  shown  in  Figure  6,  while  the  corresponding  volume 
scattering  functions  at  0.53  pm  are  plotted  in  Figure  7. 

In  the  absence  of  reliable  spectral  bottom  reflectance  measurements, 
the  reflectance  of  Florida  beach  sand,  Kansas  soil  (representing  mud), 
and  wheat  leaves  (representing  green  aquatic  vegetation)  were  used  to 
approximate  the  reflectance  of  actual  bottom  materials.  The  reflectances 
of  these  materials  in  the  range  of  0.450  to  0.575  Pm  wavelength  are  shown 
in  Figure  8 . 

Using  these  water  and  bottom  parameters,  radiances  were  calculated 
at  a^  platform  altitude  of  1 km  for  clear  sky  conditions  (23  km  visibility) 
with  a solar  zenith  angle  of  20  degrees.  Radiances  were  calculated  at 
1 meter  depth  intervals  1 rom  0 to  20  meters.  A surface  reflectance  of 
0.021  was  used,  which  is  valid  for  observation  angles  less  th<an  20°  and 
windspeeds  less  than  10  meters  per  second  l5l.  The  atmospheric  phase 
function  for  Deimend j fan' s ( 10  1 haze  M (marine  or  coastal)  size 
distribution  was  used  for  calculating  atmospheric  effects. 
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SCATTEFtING  ANGLE  (degrees) 

HC.UFtE  7.  VOLUI^  SCATTERING  FUNCTIONS  VS.  SCATTERING  ANGLE  FOR 
PETZOLD'S  STATIONS  9,  5,  AND  2200 
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' : i>l  i Ml.  I . m.1 III  (1  i s |i  1 .lyi'd  in  .st'vt'ral  w.ivs  . One 

' ' !.m.  , wlii.  h is  iisMlul  for  .inalysis  of  bottom 

■ . • ' . ii  iiiiu'l  plot  shown  in  I’ipiiros  9 aiui  10.  Tliesc 

"ri  i IS  I.. tv"  ot  t ho  siKn.iIs  .is  the  dopth  varios  from 
• i!  . : ,’  i , ;i,:  1 t.i  intinitv  (lowor  loft  end).  Figure  9 is  for 

'.".1  • . i ''  ..  t mu  |ii  i ; u.ilo.  .ittonu.ition  ooo  f f i c i en  t s , showing  curved 

ti  l l I !.■  . figuii  In  Is  tor  w.i vi‘ 1 ongi hs  with  equal  attenuation  cocf- 
ient  .;  t !■.  t t i ; i-.  t . . r i es  in  this  o.ise  are  ne.irly  straight,  with  only 
1 light  urv.ituio  duo  ti'  si.ittoring  and  internal  reflection  effects. 

-■iriil  ir  plots  I .in  ho  gonor.itod  for  actual  scanner  data,  except  that 
the  a.  tn.il  d.it  i |ioints  .iro  distributed  around  the  ideal  trajectory  due 
to  noisi  . Fipuri'  11  shows  <i  scattir  plot  of  dat.a  values  along  a transect 
in  ;>t.  Andrew  H.iv , Klorid.i,  collected  by  F.KlM's  M-7  Scanner  111).  The 
two  t I 1 jei  tor  ios  indie. itod  on  this  figure  were  calculated  for  the  water 
p.ir.imotors  measured  at  tlu‘  lime  of  the  overflight,  with  bottom  reflectance 
ratios  tor  s.ind  (lower  curve)  and  veget.ation  (upper  curve).  An  abrupt 
change  in  bottom  type  occurs  at  .about  2 meters  depth,  indicated  by  a 
shift  1 rom  one  trajectory  to  another. 

4.3  DKSCRIl’TION  OF  MRA  TF.CHN  KlUF 

The  modified  ratio  algorithm  (MRA)  uses  the  ratio  of  the  bottom- 
reflected  signals  in  two  channels  to  discriminate  between  different 
bottom  types  under  a v.iriahle  depth  of  water.  This  technique  is 
illustrated  in  Figure  12.  Curves  A and  11  in  tliis  figure  represent 
the  itlealized  "trajectories"  of  signals  for  two  distinct  bottom  types. 

A line  beginning  at  the  point  with  slope  R is  drawn  between 

these  curves:  this  line  is  the  "decision  boundary"  in  the  MRA  technique 

for  discriminating  between  bottom  tv|)es  A and  B.  The  values  L.j^  and 
are  the  "offsets"  in  channels  1 ami  ?,  which  are  normally  obtained  by 
sc.anning  over  deep  w ter.  If  a given  set  of  observed  radiances  and  L2 

fall  above  the  decision  boundary,  the  ratio  value  is  greater  than  R,  and 
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the  material  is  presumed  to  be  bottom  type  A.  All  points  below  this 
line  have  a ratio  value  less  than  R,  and  arc  presumed  to  belong  to 
bottom  type  B. 

If  the  signal  trajectories  were  strniglit  lines,  and  if  there  was 
no  noise  in  the  signal,  the  MRA  technique  would  allow  a perfect  dis- 
crimination between  bottom  types.  The  trajectories  are  nearly  straight 
if  the  water  attenuation  coefficients  in  bands  1 and  2 are  equal,  but 
there  is  some  curvature  due  to  scattering  and  internal  reflection 
effects.  Further,  system  noise  and  variation  of  the  svirface-ref lected 
signal  due  to  waves  cause  the  actual  signals  to  be  d ifstrihuted  around 
idealized  trajectories  rather  than  to  fall  exactly  on  them.  For  a given 
depth  z' , the  distribution  of  data  points  may  be  approximated  by  a 
Gaussian  distribution 


where  1.^  and  I,^  aie  the  radiances  at  the  deptii  z'  in  the  absence  of 
noise,  and  a ^ and  are  the  r.m.s.  noise  signals  in  bands  1 and  2. 

The  ellipse  drawn  in  Figure  12  represents  one  standard  deviation  of  the 
radiance  at  z = z'  due  to  noise.  The  equation  of  this  ellipse  is 


+ 


= 1 


O) 


and  the  probabili 
will  fall  inside 
ellipse  is  drawn 
as  a dotted  lino 


ty  that  a given  radiance  measurement  at  depth  z' 
this  ellipse  in  39. 3X.  Suppose  now  that  a similar 
tangent  to  the  decision  boundary  (this  ellipse  is  drawn 
in  Figure  12).  This  ellipse  Is  described  by  the  equation 
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(9) 


Using  Iho  probability  distribution  (6)  it  can  be  shown  that  tlie 
total  probability  tliat  a given  radiance  measurement  at  depth  z'  over 
bottom  type  A will  tail  below  the  decision  boundary  is 


P'  = i/2  erfc  (—  ^ 
r'2 


(10) 


where  eric  is  the  complemented  error  function.  This  is  equivalent 
to  the  probability  of  misclassif icat ion  for  tlie  MRA  with  this  decision 
boundary.  In  the  case  of  three  (or  more)  bottom  materials,  there  are 
two  (or  mt.re)  decision  boundaries.  The  probability  of  misc 1 ass i f icat ion 
for  a material  liaving  a decision  boundary  on  both  sides  is  then  equal 
to  tlie  joint  probability  of  a measurement  falling  beyond  either  decision 
boundary. 


4.4  KVAI.UATION  OF  MRA  TKCUNIQUK 

A comjniter  program  was  written  to  calculate  the  (irobabi  1 i t ies  of 
m i sc  lass i t i cat  ion  for  the  MRA  technique  with  the  bottom  materials  and 
water  types  described  in  section  3.2.  This  program  uses  the  radiance 
data  comiiuted  by  the  radiance  model  program  to  simulate  the  performance 
of  the  MRA  on  an  actual  remote  sensing  data  set.  The  r.m.s.  noise 
signals  are  computed  from  the  equation 
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wluTi.-  L is  tlu-  lotal  railiaiue  bfini;  measured  and  1,  is  tin.'  noisu- 

n 

fqu  i va  I cnl  radiaiui-  ol  the  sensor.  1‘or  tiu*  KRIM  M- 7 scanner,  the 

-2  - 1 - 1 

no  i se-i'<iu  i va  1 en  t radi.ince  is  approximately  .0(11  mW  cm  sr  .m  in  the 
hlue-^reeii  re>>ion  ot  the  spi*ctrum. 

I'hi’  jirojitram  is  designed  for  interact  iv'e  use  so  that  varioii.s  modes 
ot  operation  can  be  selected.  Among  the  user  inputs  are  the  wavelengths, 
lUisets,  .md  the  w.iter  depth  at  whicli  the  algoritlim  is  to  be  "trained". 

I'he  ratio  values  for  each  bottom  type  are  c.ilculated  at  this  depth,  and 
the  decision  boundaries  are  chosen  miciwny  between  e;n  h ad  j.icent  pair  of 
ratio  values.  Having  chosen  the  decision  boundaries  ftir  each  bottom 
materi.il,  the  "distance"  d between  the  div  ision  bound. ire  ,md  the 
trajectory  is  calculated  at  each  depth,  from  e<iu.it  ion  (9).  The  prob.ihilit 
oi  misc  1/iss  i f icat  ion  is  then  calculated  from  e<|uat  ion  ( lo)  .as  a function 
of  depth  for  each  botfiim  type. 

The  results  of  one  such  set  of  calculations  are  shown  in  Figure  li- 
I'hese  results  are  for  water  type  3,  with  w.ivelengths  selecttni  to  m.iximize 
the  separability  of  green  vegetation,  subji'Ct  to  the  constraint  that 
tile  attenuation  coefficients  must  be  t iie  same  at  both  wavelengths.  The 
results  show  that  vegetation  is  indeed  the  most  accur.ately  classitied 
bottom  type,  with  a probability  of  m i sc  1 ass  i 1 i c.it  i on  less  than  25 
pi'rcent  out  to  a depth  of  about  5 meters. 

The  relatively  poorer  performance  over  sand  and  mud  is  due  to  the 
similarity  in  the  shape  ot  their  spectral  reflectance  curves.  This  can 
be  seen  in  Figure  10,  which  shows  th;it  the  signal  trajectories  for  sand 
and  mud  .are  very  close  together  for  these  wavelengths.  This  limit.ation 
of  the  MRA  technique  will  be  further  I'iscussed  in  the  following  section. 
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Tho  performance  illustrated  in  Figure  1!  is  obtained  if  the  water 
quality  is  uniform  throughout  the  scene.  In  actuality,  water  quality  is 
likely  to  he  non-uniform  in  the  coastal  environment  because  of  wave 
ait  ion,  runoff  from  tlie  land,  and  biological  activity  nearshore.  In  the 
typic.il  situation  tiie  water  quality  is  similar  to  water  type  1 in  deep 
water,  and  to  water  type  3 or  even  water  type  5 in  shallow  water 
(although  water  type  5 is  usually  more  typical  of  inland  waters).  The 
effect  of  these  ch.inges  in  water  quality  is  two-fold.  First,  since 
the  oltsets  are  usually  obtained  by  scanning  over  deep  water,  these 
offsets  will  be  incorrc'ct  if  tlie  water  quality  in  deep  water  is 
different  from  that  in  shallow  water.  Second,  since  the  optimum  wave- 
lengths and  decision  boundaries  vary  with  water  quality,  an  algorithm 
which  has  been  optimized  for  one  water  type  will  perform  less  well  in 
another  water  type. 

The  effects  of  such  changes  in  water  quality  are  illustrated  in 
Figures  14  and  15.  Fi gure  1 b shows  the  MRA  performance  for  vegetation 
in  three  situations.  Curve  A is  for  uniform  (type  3)  water  quality 
throughout  the  scene  (i.e.,  the  same  situation  as  in  Figure  13.  The 
offsets  wi'te  obtained  from  water  type  3 and  tho  operating  parameters 
(wavelengths  and  decision  boundaries)  were  optimized  for  the  same 
water  type.  Curve  B shows  the  performance  in  water  type  3 with  offsets 
obtained  from  water  type  1,  but  with  the  algorithm  otherwise  optimized 
for  water  type  3.  Curve  C shows  the  results  of  the  same  algorithm 
applied  to  water  type  5:  that  is,  the  offsets  were  obtained  from 

water  type  1,  the  algorithm  was  "trained"  in  water  tvpe  3,  and  the 
performance  was  evaluated  in  water  type  5.  Figure  15  shows  the 
average  performance  for  the  three  bottom  types  in  the  same  three 
situations. 

The  performance  in  shallow  water  is  approximately  as  expected. 

The  average  probability  of  ml  sc  lass i f Icat ion  is  the  smallest  for  the 
optimized  algorithm  (curve  A)  out  to  about  7 meters  depth.  The 
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FIGURE  13.  PROBABILITY  OF  MISC LASSIFICATION  VS.  DEPTH  FOR  WATER  TYPE  .5 
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per t ormance  in  situation  h is  lortuit iously  somewiiat  better  for  vegeta- 
tion, but  is  worse  lor  sand  ami  mud.  The  perfornuince  in  situation  C is 
also  worse  than  the  optimum  case,  as  expected.  Somewhat  surprisingly, 
liowever,  the  average  proh.ibility  of  misclassif  ication  is  rehitively 
high  (17  ) even  at  zer('  depth.  This  suggests  a fundamental  limitation 
ol  the  MKA:  namely,  that  the  bottom  discrimination  technique  used  in 

this  algorithm  is  not  an  optim.il  one,  since  it  incurs  a relatively 
l.irge  error  even  in  the  absence  of  an  overlying  water  column.  There 
are  two  reasons  even  in  the  absence  of  an  overlying  water  column. 

There  are  two  re.isons  for  this  error.  First,  the  choice  of  wavelengths 
is  constrained  by  the  requirement  of  equal  .attenuation  coefficients, 

.and  the  wavelengths  cannot  therefore  be  optim.ally  selected  for  bottom 

[ discrimination.  Second,  since  the  metliod  uses  ratios  of  signals,  bottom 

r 

types  with  nearly  equal  reflectance  ratios  (such  .as  s.and  and  mud)  are 
not  well  separated.  The  development  of  algorithms  not  subject  to  these 
limitations  is  discussed  in  section  3.5. 

The  performance  of  the  MRA  in  deep  water  is  perhaps  counter- 
intuitive, but  can  be  rpiite  easily  explained.  For  the  optimized 
algorithm,  the  proh.ibility  of  misclass i f icat ion  approaches  100  percent 
I for  .all  bottom  types  as  the  depth  .approaches  infinity.  This  is  to  be 

expected,  since  the  bottom  cannot  he  recognized  correctly  If  it  cannot 
be  seen.  Wlien  the  proper  offsets  are  used,  the  numerator  and  denominator 
of  the  ratio  both  approach  zero,  and  the  ratio  value  is  dominated  by 

I 

noise:  i.e.,  the  value  becomes  purely  random.  Wlien  an  Improper  offset 

is  used  in  the  denominator  (as  in  cases  B and  C) , the  ratio  value 
approaches  a finite  limit  as  the  depth  Increases.  Therefore,  the 
probability  that  the  bottom  will  be  classified  into  one  or  more  of  the 
; categories  approaches  a non-zero  limit,  even  though  the  bottom  cannot 

Ibe  seen.  The  performance  of  the  MRA  thus  appears  to  be  better  for  the 
non-optimum  cases  than  for  the  optimized  algorithm  In  deep  water. 

In  practice,  probabilities  of  ml  sc  lass  if icat ion  greater  than  67  percent 
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(tor  a throe  catoKorv  situation)  are  not  of  interest  anyway,  since 
this  is  tile  ()rohali  i 1 i ty  ol  m i se  1 ass  i f icat  ion  for  a purely  ramioni  guess. 

A summary  of  the  MRj\  performance  in  tlie  t It  rei'  water  types  is 
shown  in  Table  3.  Tlie  number  in  the  first  column  indii-ates  the  water 
tvpe  in  whielt  the  evaluation  was  performed.  The  second  column  indicates 
the  water  type  from  whicli  tlie  offsets  were  taken,  and  the  third 
c-olumn  indicates  the  water  type  for  which  the  algorithm  was  optimized. 

I'he  last  two  columns  contain  the  average  probability  of  m isc  1 ass  i f i cat  ion 
at  zero  depth,  and  the  depth  at  which  this  probability  is  equal  to 
50  percent,  respectively. 

The  first  three  entires  in  Table  3 represent  fully  optimized 
cases,  in  the  sense  that  offsets  were  chosen  from  the  same  w.ater  type 
and  wavelengths  were  chosen  to  equalize  the  attenuation  coefficients. 
However,  the  wavelengths  used  for  water  type  1 were  apparently  more 
optimum  than  those  for  water  type  3 from  the  point  of  view  of  bottom 
discrimination,  since  the  probability  of  misclass i f icat ion  at  zero 
depth  is  smaller  for  type  I than  fc^r  type  3.  Similarly,  the  wavelengths 
for  type  3 are  better  than  those  for  type  5 for  bottom  discrimination. 
This  explains  why  the  MRA  performance  in  water  type  5 using  the  type  3 
.algorithm  (the  last  case)  is  better  than  the  performance  using  the 
"lully  optimized”  type  5 algorithm. 

TABl.K  3 

SUMMc\RY  OF  MR.\  I’KRF ORMANCF. 


WAiKR  TYJ’i: 

OFKSFiS 

AI.CORITIIM 

l’m(^^0) 

z(Pm  = ' 

1 

1 

1 

1 3 . 2% 

K . A m 

3 

3 

3 

lb.  6^ 

6.0  m 

5 

5 

5 

26.7% 

3.0  m 

3 

! 

3 

17.2% 

6.0  m 

5 

1 

3 

17.2% 

3.6  m 
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4.5  KURTHKR  AI.CUR  ITll.M  DKVKI.OPMl'.NT 

Tlif  prim.,  rv  limit.!  Lion  ot  tlie  MRy\  torhniqiio  is  its  in.ibi  1 i ty  to 
sop.ir.ito  bottom  m.itiri.ils  with  ecjua  i or  lUMrly  fqu.il  v.ilufs  of  the 
i)ottom  rft’lfit.iiu'f  r.itio  in  tlio  two  oper.iting  channels.  For  such 
bottom  tvpes,  the  siRn.il  t r.a  jec  L or  i es  are  coincident,  or  nearly  coinci- 
dent, .IS  in  Fi>;ure  10.  This  results  in  the  relatively  bi^h  probabilities 
('I  mi  sc  1 .iss  i I i I .1 1 ion  lor  s.ind  .ind  mud  shown  in  Figure  13. 

It  is  possible  to  select  wavelength  pairs  with  distinct  sign.al 
L ra  ji-ctories  lor  any  set  of  bottom  types,  if  the  attenuation  coefficients 
are  not  required  to  be  equ.il.  i he  trajectories  in  Figure  9,  for  example, 
are  distinct  evim  though  the  ratios  of  bottom  reflectances  are  nearly 
the  sami'  lor  s.ind  and  mud.  In  this  case,  however,  the  trajectories 
c.inni't  be  sep.ir.ited  by  ,i  1 ine.ir  decision  boundary  as  in  the  MRA. 

An  .ilternative  processing  .algorithm  for  this  case  is  suggested 
by  the  .ijiprox  i m.i  te  1 y I'xponential  relationship  between  signal  levels 
.ind  w.iter  di-pth.  If  Llie  transformation: 

V.  = In  (I..  - I,.  ) 

1 1 IS 


is  applied  Li-  the  d.it.i  .and  the  results  .are  plotted,  the  result  is  a set 
ol  ne.arly  linear,  par.allel  t r.a  jec  t or  i es  . Figure  lb  shows  such  .a  plot 
tor  the  same  d.ata  set  .as  plotted  in  Figure  . The  slope  of  these 
trajectories  is  given  by  the  r.atio  of  the  .attenuation  coefficients  In 
the-  two  ch.annels.  Thus,  .a  linear  decision  bound.ary  can  be  dr.awn  between 
each  pair  of  trajectories  to  separate  the  bottom  types.  These  decision 


bound.aries  .are  described  by  the  equ.ation 
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Itu’  "d  i s t anil"  bfiwi'fii  tlie  t ra  joe  tor  i t-s  ol  material  A and  material 
B are  Kiveii  apiprox ima te  1 y by  the  equation 


S 


Ki  + 


(14) 


where  is  the  bottom  refloctanee  for  material  A in  band  1,  etc. 

ibis  equation  is  approximate  because  it  neglects  scattering  effects. 
However,  such  an  index  allows  a criterion  for  optimum  band  selection 
for  separating  any  given  set  of  bottom  materials. 

This  algorithm  holds  promise  because  it  is  not  restricted  to 
wavelength  pairs  with  ecjual  attenuation  coefficients.  Thus,  the 
wavelengths  can  be  chosen  for  optimum  discrimination  of  the  bottom 
types  in  the  siene.  Further,  the  method  is  not  inherently  limited  to 
two  channels,  .as  the  MRA  is.  It  is  therefore  expected  that  better 
pertorm.ince  could  be  obtained  with  this  method  than  with  the  MRA, 
especi.illy  in  the  cast.’  of  materials  with  similar  bottom  reflectance 
ratios.  We  hope  to  carry  out  an  error  analysis  for  this  algorithm 
during  the  next  year,  as  well  as  to  include  scattering  effects  in 
equation  (11)  and  define  a criterion  for  optimum  channel  selection. 
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A.  I DESCRIPTION  OF  PROGRyVMS 
A. 1.1  ABSCAT 

The  ABSCAT  program  calculates  the  absorption  and  scattering 
coetficients  tor  a mineral  wlien  given  the  experimental  t ransmi t t ances 
and  reflectances.  The  user  invokes  this  program  by  typing 

$RUN  ABSCAT  1 = CONTROL  FILE 

where  the  control  file  contains  a namelist  input.  An  example  of  this 
control  file  follows: 

&LIST 

RAP  = .6,  RPA  = .03.  PU  = 5*1.,  PL  = 5*0,  DP  = 5*. 01,  SIC  = , 
ITER  = 20,  NPAR  = 5, 

6.  END 


These 

parameters 

are  defined  as: 

IO\P 

- related  to 

index  of  refraction 

RPA 

- related  to 

index  of  refraction 

PU 

- parameter 

upper  bounds  (AL,  FS , 

S) 

PL 

- parameter 

lower  bounds  (AL,  FS , 

S) 

DP 

- parameter 

increment  for  curve  fitting 

SIC 

- parameter 

error  allowance 

ITER 

- number  of 

iterations  for  Newton 

's  method 

NPAR 

- number  of 

par.imeter  bounds  used 

The  program  will  first  prompt  the  user  to  declare  mineral  input  and 
coefficient  output  files.  The  next  prompt  will  be  for  the  user  to 
enter  the  sample  number  and  thickness  for  the  infinite  thickness  sample. 
If  there  is  no  infinite  thickness,  the  user  enters  zero  for  the  sample 
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numbiT  .ind  lla-  lliickiu'ss  of  tlie  tirsl  s.implc.  Tlie  next  pri;mpL  is  for  tip 
to  ttirii-  t ransm  i 1 1 ani’o  sample  numbers.  The  iisiT  will  t lien  be  prompted 
lor  till'  eor respond  1 ng  reflectance  numbers.  The  final  prompt  will  be 
for  till  correspond  ing  sample  t b i r knosses . Note  that  11  there  was  not 
any  infinite  thickness  sample,  the  first  sample  thickness  was  already 
entered.  Therefore,  the  remaining  thicknesses  shoulii  be  entered. 

A.  1.2  I'KKDirr 

The  I’KKDICT  program  calculates  predicted  transmittance  and  reflect- 
ances given  the  scattering  coefficients  and  the  mineral  thickness  for 
which  the  prediction  is  to  be  made.  ihe  user  invokes  this  program  by 
typing 

$RU.N  I’KKDICT 

The  user  will  then  be  prompted  to  enter  the  output  file  and  the  input 
file  which  contains  the  scattering  and  absorption  coefticients  (output 
from  AB.SCAT).  The  next  |>rompt  will  be  lor  the  mineral  thickness  at 
which  to  make  the  prediction  and  RAP  and  RPA. 

A.l.T  OPAgilK 

The  Ol’AQUK  program  will  take  experimental  ref  lect  anci-s  lor  a 
mineral  which  were  determined  foi  different  thicknesses  and  smooth  this 
data  using  equation  (3).  The  average  reflectance  over  all  thicknesses 
will  be  found.  The  program  will  then  output  a file  compatible  with  the 
output  of  PKKDICT  that  contains  the  smoothed  and  averaged  ref  K'Ct  .inces 
and  all  t ransml t t ances  set  to  zero.  The  user  will  be  prompted  tor  tlu' 
input  file,  output  file,  and  experimental  t h i cktusses . I'p  1 1'  four 
thicknesses  may  be  used.  The  user  invokes  this  program  bv  tvpmg 


$RUN  OPAQUK 


'erjm 
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A. 1.4  STAIN 

riu'  prcigrum  STAIN  will  t.ikf  an  input  tile  containing  one  .set  of 
exper  imi'iit  a 1 t ransni  i t t.ince  and  reflectance  spectra  and  smooth  the  values 
ai'cording  to  equation  (3).  The  output  will  be  compatible  with  the  out- 
put of  till'  I'RKDICT  program.  'The  user  invokes  this  program  by  typing 

$RIIN  STAIN 

A.  1.5  SANDRKF 

The  SANURKK  program  will  model  sand  given  the  predicted  transmit- 
tance and  reflectance  files  (i.e.,  outputs  from  I’RKDICT,  OPAQUK,  and 
STAIN)  of  the  constituent  minerals,  the  sand  configuration  (e.g.,  particle 
size),  and  the  sensing  configuration.  The  user  invokes  this  program  by 
typing 


$RUN  SANDREK 

The  user  will  be  prompted  for  an  output  file  which  will  contain  the 
I predicted  rellectance  spectrum.  The  user  will  then  be  prompted  for  an 

input  file.  The  input  file  will  contain  the  names  of  the  files  contain- 
ing the  predicted  transmittance's  and  reflectances  which  will  make  up 
the  sand.  The  file  will  also  contain  a series  of  namelists  which  will 
describe  various  experimental  configurations  of  sand  composition  and 
illumination-sensor  geomet ries . Table  A. 1 defines  these  parameters  and 
T.ible  A.  2 gives  an  example  of  the  input  file.  The  data  for  Quartz  is 

defined  as  90.2%  transmittance  and  9.8%  reflectance  over  all  wavelengths. 
Therefore  if  Quartz  is  listed  as  a mineral  file,  a file  need  not  exist 
and  the  program  will  assign  tlie  above  values.  The  user  will  also  be 
prompted  as  to  plotting  tlie  reflectance  spectrum  for  the  sand. 

1 


60 


I lUUf 


t ' 'A 


A fti  Utk  1 I )Hit 


tut  ■ f MS'T  » Of  MU 


• AN 


ERIM 


r.ihlo  A.  1 . SANDRl.K  IniiuL  I’araini-t er  Dot  ini  t ions 

N’MIN':  N'umhor  ol  minora  Is  in  saiui 

N(;iA)M:  N'umbor  oi  goomolry  namolisls 

X(l):  DoptI)  of  lop  layor  of  sand  (negative  number) 

X(d):  Itoplli  of  b.)ltom  layi'r  of  sand  (negative  number) 

I’H  ! ( ! ) -I’ll  I (i>)  : Polar  view  angles 

Til:  Sun  (illuminator)  angle 

l’SI(l)-l’SI  7):  Asimuth  view  angles 

NAZ : Number  of  azimuth  view  angles 

Nl’Ul.:  Number  o!  polar  view  angles 

Sf 1 , 1 , 1) -S(2 , 2 ,b) ; HNV  values  for  minerals 

deli  nit  ion:  H = '■ ' n V = 1.414*11 

the  array  storag-  lor  S(l,  .1,  K)  is  defined  by 


.1  1 

J2 

L 

9 

I 

2 

I - 

2 - 

11 

V 

H 

V 

r-l  , 

r:  M * 

O 5 - 

.':l  ^ " 

To|) 

Bof  t-om 

Layer 

Layer 

NSKY:  = 0 if  no  sky  (direct  sun) 

=1  if  sky  (diffuse  illumination) 

I ■ 
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r.ibU-  A. 2.  Kx.imp  li'  ot*  SANDRKF  Input  Kilo 


I I,  (NMIN) 

n I »,  1 . (MINERAL  FILE) 

( Mt  1 A.  t-n.  ( I , in  (BACKGROUND  FILE) 
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